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We study the cosmological perturbations of the new bi-metric gravity proposed by Hassan and 
Rosen [1] as a representation of massive gravity. The mass term in the model, in addition of ensuring 
ghost freedom for both metrics, causes the two scale factors to mix at the cosmological level and 
this affects the cosmological perturbation of the model. We find two combinations corresponding 
to the entropy and adiabatic perturbations of the theory. In this sense we show that the adiabatic 
perturbations could be a source for the entropy perturbations. So in addition to the adiabatic 
perturbations, entropy perturbations can also be present in this theory. We also show that the 
adiabatic perturbations are not constant at the super horizon scales, implying that the theory could 
not be used to describe the inflationary epoch, even if it can impose some corrections to the standard 
inflationary scenarios. 
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I. INTRODUCTION 



The current problems facing standard cosmology are 
of such breath and depth that any fixation to these prob- 
lems requires a deeper understanding of the prevailing 
theory that has been and is still in use, namely the gen- 
eral theory of relativity. Novel ideas, both old and new, 
have been suggested over its rather long history in the 
hope of a partial remedy to some of the most pressing 
problems we are facing today. Two interesting observa- 
tions made over the period of the past two decades, that 
is, the accelerated expansion of the universe and galaxy 
rotation curves, amongst others, have stirred a plethora 
of activities aimed at modifying the standard Einstein- 
Hilbert (EH) action and hence to offer a solution to some 
of these problems. With partial success, one may con- 
sider various modifications to the standard EH action in 
the form of modified theories of gravity with and without 
torsion, gravity in extra dimensions, and more recently, 
massive gravity, to name but a few. The latter has been 
attracting the attention of many experts in the field in 
the past few years. It is based on an old idea by Fierz and 
Pauli (FP) [2] where an effective field theory with a mas- 
sive graviton was proposed. This theory was, however, 
problematic at the linearized level since the correspond- 
ing Newtonian potential was discontinuous for a vanish- 
ing mass, m? , resulting, for example, in a large correction 
to the deflection of light around the Sun compared to 
the experimentally accepted value [3] predicted by Gen- 
eral Relativity (GR). This is referred to as the vDVZ 
discontinuity. The source of the discrepancy was traced 
to the degrees of freedom of the graviton, being two for 
a massless and five for a massive graviton. Another set- 
back was discovered later on when it was demonstrated 
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that when self interacting terms are added to the action, 
ghosts would appear in the theory [4]. It has only been 
in the past few years that a method for fixing the above 
problems has been realized [5]. Theories rooted in the 
FP action are collectively known as massive gravity. 

Perturbation methods are and have been an integral 
part of attempts to find solutions to complicated prob- 
lems. Massive gravity is therefore no exception in this 
regard. In non-linear theories such as GR, perturbation 
methods can be particularly useful when one seeks the ef- 
fects of small changes in the metric. Such methods have 
been exploited in recent years to study, for example, the 
structure formations in the universe. It then seems only 
natural to conduct such a study when dealing with mas- 
sive gravity. Such a study becomes even more attractive 
if one considers the fact that massive gravity is inherently 
a bi-metric theory. This comes about since any modifi- 
cation to the EH action in the form of a self-coupling 
term involving no derivative whose definition is based on 
one particular metric requires an additional metric which 
may be dynamical or fixed [6, 7]. The appearance of a 
second metric in the theory can be understood on general 
grounds. If the second metric is non-singular, spherically 
symmetric, and both are diagonal in some coordinate sys- 
tem, a Killing horizon for one metric must be a Killing 
horizon for the other [8] . Interestingly, it has been shown 
that the off diagonal elements undergo no modification at 
large distances [9] . From a cosmological perspective, in a 
bi-metric massive gravity theory with the second metric 
static, there is no spatially flat FRW solutions [10], con- 
trary to the bi-gravity formulation of the FRW cosmology 
for which homogeneous solutions exist [11, 12]. 

It has increasingly been realized that modern obser- 
vational date can be explained by perturbation methods 
which have become an integral part of any study deal- 
ing with cosmological fluctuations. On the other hand, 
from a theoretical point of view, the inflationary scenario 
can also be supported by the modern observational data 
rather accurately. This makes the studies of the cosmo- 
logical perturbations in massive gravity [13] all the more 
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interesting since its predictions are gradually becoming 
known and are still not on firm grounds. As was men- 
tioned above, one realization of massive gravity is in the 
language of bi-metric theories. Consequently one may be 
interested in a cosmological perturbation theory in such 
models, for the existence of two metrics in these theories 
may result in non-trivial features, specially in the behav- 
ior of relative metric perturbations. In this paper we will 
consider scalar perturbations for both metrics in a mas- 
sive gravity model. We will obtain two gauge invariant 
combinations of the perturbed functions which can be 
responsible for the adiabatic and entropy perturbations. 
As a result we will obtain the super horizon limit of the 
equations of perturbations and show that the adiabatic 
perturbation is not constant in this scale. 

The scope of the paper is as follows: in the next sec- 
tion we present the bi-metric model studied here. Section 
III deals with the the cosmological perturbations of the 
model and definition of the gauge invariant quantities. In 
section IV we obtain the equations of motion for the adi- 
abatic perturbations and discuss the super horizon limit. 
Conclusions are drawn in section V. 



II. THE MODEL 



in the metrics / and g in the sense that 



We begin with the bimetric action [1] 
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where R{g) and R{f) arc the Ricci scalars correspond to 
metrics p^i/ and /^i, respectively, and Pi are some arbi- 
trary constants. We define the polynomials e„ as 
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Without the kinetic term for metric the action de- 
scribes a theory of massive gravity which is ghost free to 
all orders in the decoupling limit [5]. The cosmology of 
such a theory was studied in [10]. 

The main point of the mass term is that it is symmetric 



/ d^x^gY^fi^e^(y^) 

n. 4 

= / d4a;y^^/3„e4-„ (/Pg) ■ (4) 



So one can read off the equation of motion for /^j^ by 
having the equation of motion for g^^i,. The equation of 
motion for g^j^y can be easily obtained, resulting in [14] 
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where is the Einstein tensor for metric g^j^. 
The equation of motion for /^^ is 



0, 
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+ a2 {Yl, - Vei(Y) + ^.e2(Y)) 
- a3 {Jl, ~ Y^,ei{Y) + Y^,e2{Y) - ^.e3(Y)) , 

(8) 

where F^i, is the Einstein tensor for metric /^^ and we 
have defined Y = f~^g and = /34_„. The cosmolog- 
ical solutions of this bi-metric theory has been considered 
in [11]. 

III. COSMOLOGICAL PERTURBATIONS 

In order to study the cosmological perturbations of the 
model, we consider the scalar metric perturbations of the 
two metrics (our notation is compatible with that in [15]) 

dsl = al{t)[- (1 4- 2(j)i)dt^ + 2d,Bidx'dt 

+ [(1 - 2T/'i)(5y- + 2d^djEi]dx'dx^] , (9) 

and 

ds] = al{t) [ - (1 + 202)dt^ + 2diB2dx'dt 

+ [(1 - 2tlj2)kj + 2d,djE2Wdx'] . (10) 
Let us define four zero order quantities for metrics g^j^u 

Xi= l3iala2 + 2P2aial + 13-ial, (11a) 

A2 2/3oaJ + 7pia\a2 + 8l32aial + Zp^al, (lib) 

A3 = /?oaJ + 2/3ia^a2 + /?2aia^, (lie) 

A4 = Poal + 3/3ia^a2 + 3/32aia^ + P^al- (Ud) 
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Y2 = 2l3i4 + Ipzalai + 8/^2020? + Wial, (12b) 
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These terms can simplify the mass parts of the equations 
of motion. The background equations for metric g^j^i, is 
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where ' = d/dt and Hi — a'l/ai. Similarly for f^^ we 
have 
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where i72 = 02/02- For the Kinetic term of the metric 
g^ui we obtain, to first order [20] 

G^^ ^ 2V2 fv'i + ifi(£^i - Bi)] - 6Hi^[, (15a) 
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+ eV'i + v^D, (i5d) 



with 
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The mass term then becomes, to the first order (sec the 
Appendix) 
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The kinetic term for the metric f^^, is represented by its 
corresponding Einstein tensor F^^ which can be deduced 
for the first order perturbations by making transforma- 
tion 1 — > 2 in equations (15a)-(15d). The mass term 
can be obtained easily by transformations 1 -f-?- 2 and 
/5ri — > Pi^n ill equations (17a)-(17d). The result is 
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As we know, general relativity is invariant under coor- 
dinate transformations. Hence we restrict our model to 
scalar perturbations so we only consider scalar coordinate 
transformations as follow 



t^t + St, 
x'' ^ + 6'-'djSx, 



(19a) 
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and consequently the scalar metric perturbations behave 
like [21] 
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where i = 1,2. We then define six independent gauge 
invariant quantities as follows 
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with i = 1,2. Using (21c) and (21d), one can define a 
gauge invariant quantity 



A = £' -B, 



(22) 



(26b) 



4l^-^^5.(f'-A), 
4' = ^^^9.9,8, (26c) 



which simphfies the calculations that follow. In order 
to write the field equations in the gauge invariant form 
we can fix two scalar gauge freedoms by the following 
conditions 



Ei+E2 = 0, Bi+B2^ 0. 



(23) 



This gauge fixing has the advantage that the symme- 
try between the two metrics remains true, even after the 
gauge fixing. We then rewrite the kinetic term to first 
order for metric g^^i, i.e. relations (15a)-(15d) in terms of 
six gauge invariant quantities ^'1,2, 'i'1,2, A and £ as 
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(24b) 
(24c) 
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6(i72 + 2i/()($i + *i), (24d) 



and similarly for metric /^i^ we obtain [22] 
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= [2(^2 - ^2) + [Hi + 2H'^)£] 



i=l 



+ a*^' + 6i72($2 + 2^-2)' + 'i{H'^ + i?2i?2)A 
+ 6(i72 + 2i/^)($2 + *2). (25d) 

In terms of gauge invariant quantities, the mass terms 
(17a)-(17d) can be rewritten as 
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and similarly for (18a)-(18d) 
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By transforming to a spatial Fourier space, is replaced 
by —k'^, where k is the wave number. Using the above 
equations and after some algebra we arrive at writing £ 
and A as functions of 4>i and and deduce the following 
constraints 
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J_ _ ^1 + X2 ^ /r^2) The second constraint can be simplified with the use of 

Ci " Aai 2 ^' (28) and (35) with the result 



and 
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Using the background equations (13a), (13b), (14a) and 
(14b), one can simplify the quantities Ai and Ci as follows 



Ai = ]-Ci = '^{2Xi-X,-X2)-\ (35) 
2 nig 
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^ "^Z namical equations then become 



A2 = -C2 = ^^{2Y4 — Yi — Y-))^^ , (36) which can be used to simplify the calculations. The dy- 



~2kHi - 3i7i(2*i - HiA)' 



('2(*2 - *!) + (i/i + H2)A - 



ai 



('2$i +i?iA + A') 



(38) 



-2kH2 - 3/^2(2^-2 + i?2A)' = TO/ 
and 



- -^Yi f2(«'i - *2) - {Hi + H2)A + h^s) + — (2$2 - -ff2A - A') 
2a2 V 3 / 02 ' 



[2($i - *i) - {Hf + 2H[)£] + 6Fi($i + 2*1)' + 6*'/ ~ 3{H'^ + HiH[)A + 6{H^ + 2H[){<i>i + *i) 



ai 



(39) 



k^iX-i - Xi)E - 6(X3*i + Xi^2) + KHiXz - H2Xi)K + -Xi [2($2 - *i) - [Hi + i?2)A - 2A'] 



(40) 



- fc2 [2($2 - *2) + [Hi + 2H'2)£] + 6i?2(*2 + 2*2)' + &-^2 + 3(^2 + H2H'2)K + &{Hl + 2F^)($2 + *2) 



02 



k^{Y:i - Yi)8 - 6(^3 ^-2 + n^i) - 3(i72i^3 - HiYi)K + -Yi [2($i - $2) + {Hi + i/2)A + 2A'] 



r 



(41) 



where we keep in mind that the quantities A and £ must 
be replaced by (28) and (29). In this sense we obtain 
four equations for four dynamical variables. 



IV. EQUATIONS OF MOTION AND THE 
SUPER HORIZON LIMIT 



From [12] wc know that a classical solution for the scale 
factors can be deduced by assuming oi =02. Due to this 
special solution one may expect that a suitable combina- 
tion of the variables ^ = '^i + "^2 and $ = $1 + $2 is 
responsible for the adiabatic perturbations. Similarly, a 
suitable combination of the variables ij} = i — ^^2 and 
(/) = $1 — $2 which is orthogonal to the adiabatic one is 
responsible for the entropy perturbations. To this end. 



we assume oi = 02 = a. We also define constants Xi and 
Vi as 

X,=a^x,, Y, = a^y„ (42) 

which, using equations (13a) and (14a), imply that 

mgX4 = mfUi. (43) 

In order to write the equations in terms of the new vari- 
ables ^E", (f>, and (j), we define 
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The constraint equations (28) and (37) can then be writ- 
ten as 



a(* - $) + /3(V' - (/)) = 0, 
a(i?$ + *') + /3{H(j} + V/) = 0. 

Equations (38) and (39) take the form 
13' 
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Also one can rewrite the relation (29) as 
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Equation (51) has an important property, showing that 
the entropy perturbation can be a source for the adiabatic 
perturbation. 

At this stage it is appropriate to take the following 
steps to make the equations less complicated and there- 
fore easier to solve. From the first condition (46) we have 



(53) 



and by the weighted addition of equations in (51) and 
(52) one can get 



{a^ + /3V)' + (a$ + /30)' _ (a$ -)- /30) (54) 



2H- — ]{a^ + /3V0 = 0. 



(55) 



Now by plugging (53) into the above equation we have 



- ^ ) (a* + Pi') = 0. (56) 



The second condition (47) results in 

(a*-f/3V)' = -i?(a* + /3(/)), (57) 
and reduces to 

(a«' + /3?/')' + -ff(a* + /3V) = 0, (58) 
due to (53). Comparison of (56) with (58) shows 
fc2 



3H 



(a* + ;3V) = 0, 



which leads to 

a^f = -/Sip, 
and consequently from (53) 

acf) = -d(l). 



(59) 



(60) 



(61) 



We note that (60) and (61) render equations (51) and 
(52) equivalent. So by plugging the above relations into 
(51) one finds 

(«- + $)'-—(* + $) + 27i71' = 0, (62) 



with 



7 = 1- 



3a 



a2 -/32 



(63) 



The same procedure as above can be employed for the 
remaining equations (40) and (41). They, in conjunction 
with equations (60) and (61), can be used to arrive at an 
equation of the form 

+ $" + H{^' + $') + 2-/H'{2^ - $) = 0, (64) 

which is independent of (62). So for four variables ip, 
(j), and <i> we have four independent equations (60), 
(61), (62) and (64). The two latter equations can be 
converted to a first order differential equation for and 
<f> by differentiating equation (62) and subtract the result 
from (64). The result is 

^ + ifV* + $)' -f'^ + 2-iH'\ (* + $) 



2-iH^' + A'^HH = 0. 



(65) 



Equations (62) and (65) constitute our final resulting 
equations for ^ and $. Now, combining the above equa- 
tion with (62) results in 

{^\, + $)" -2H{^ + $)' + (^y - 27i?2^ (^i + $) = 0, 

(66) 

which is the same as what is deduced in [17]. Since H = 
— 1/t as a consequence of the background equations, the 
solution of the above equation is 



f + $ = Cijn 



kt 



kt 



(67) 
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where ci and C2 are integration constants, jn{x) and 
yn{x) are spherical Bessel functions and 

^= 2 ■ 

Consequently, using the above solutions and either of 
(62) or (65), one can find '^{t) and $(i) separately. This 
means that by employing (60) and (61) one can find il){t) 
and (j){t) and eventually ^i, ■02, 0i and 02 as functions 
of time. 

For super-horizon modes, i.e. ^ 0, equation (66) 
results in 

* + $ = Cl^ + C2^"^ (68) 

and by considering the fact that t is the conformal time 
and t G {— oo,0}, the term "cit" is the damping term. 
The growing mode, i.e. "02^"^" term, results in 

ii, = -C2-t-^, q> = c2^-^t-^ (69) 

7 7 

This result suggests that in the theory of bi-metric grav- 
ity we have studied in this paper, the curvature pertur- 
bation is not constant at the super horizon limit which is 
compatible with [17]. The mass term in the theory can- 
not then be used as a manifestation of inflation. How- 
ever, this result is acceptable because massive gravity is 
responsible for the late time accelerated expansion of the 
universe. So, if we add some sort of inflationary scenario 
to the theory, one may expect the above result to impose 
some corrections to the curvature perturbations at the 
super horizon limit which restricts the parameters of the 
theory to that of the slow roll parameters. 

V. CONCLUSIONS 

Modern observational data have become so accurate 
that meaningful comparisons with theoretical predictions 
are now a reality. This has led to the introduction of var- 
ious calculational techniques, including metric perturba- 
tion methods, which have been playing an increasingly 
important role. In this paper we have employed such 
a method to study a bi-metric massive gravity theory. 
The present work proposes a full discussion of the scalar 
perturbations of the bi-metric gravity, with all the poten- 
tial terms at our disposal. Because the mass term mixes 
the scale factors of the two metrics, we expect both the 
entropy and adiabatic perturbations to be non-zero and 
depend upon each other. This can also be seen from the 
dynamical equations of the model. As a result, similar to 
the double scalar field models [18], one may expect that 
the adiabatic perturbation should be along the path of 
the classical solutions in the phase-space and the entropy 
perturbation be orthogonal to it. Due to [12] a classical 
solution in the framework used here is ai 02. So it 
seems it is natural to think that sum of the two scalar 
perturbations corresponds to adiabatic perturbation and 



the entropy perturbation corresponds to the difference of 
the two scalar fields. Since we have four scalar perturba- 
tions, the suitable combinations are the terms '$1 ± 
and $1 ± $2- From equation (51), one can see explicitly 
that the entropy perturbations are a source of adiabatic 
perturbations which is in agreement with double scalar 
field models [18]. 

Interestingly, equation (68) shows that the super hori- 
zon adiabatic perturbations are not constant for the bi- 
metric model discussed in this paper. We may expect 
that the mass term for graviton can only impose some 
corrections to the results of the inflationary scenario 
which we must add to the theory using other methods. In 
this scenario, the mass term in the model gets restricted 
by the slow roll parameters. This can be the subject of 
a future work. 

Note added 

During the completion of the present work, a study of 
the same subject appeared in [19]. 
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VI. APPENDIX 

In this appendix we derive the first order approxima- 
tion to matrix X = a/ g~^,f. Using equations (9) and (10) 
we obtain, to the first order in perturbation, for matrix 

G = G("'+eG'^)+C'(e^), (70) 

where 

(«'°')>fe)'«- <") 

and 

(g«)°^-2(^)\0,-02), (72) 

(^^)%,(i3i-B2), (73) 

(g^'^)^. = -2 (^^) [(Vi - ^2)S] + d,d,{E^ - E2)]. 

(74) 

We finally expand the matrix X = \/G- 

X = X("' +eX(^) +0(6^), (75) 

where 

(X(°))'^ = ^5i;, (76) 

V ' V fli 
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and 



(77) 



The expansion of matrix Y = ^ f^^g can be easily ob- 
tained by observing that metrics g^j^ and /^^ can be 
transformed into each other by transformation 1 — > 2. 
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